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We construct Holstein-Primakoff Hamiltonian for magnons in arbitrary slowly varying spin 
background, for a microscopic spin Hamiltonian consisting of ferromagnetic spin exchange, 
Dzyaloshinskii-Moriya exchange, and the Zeeman term. The Gross-Pitaevskii-type equation for 
magnon dynamics contains several background gauge fields pertaining to local spin chirality, in¬ 
homogeneous potential, and anomalous scattering that violates the boson number conservation. 
Non-trivial corrections to previous formulas derived in the literature are given. Subsequent map¬ 
ping to hydrodynamic fields yields the continuity equation and the Euler equation of the magnon 
fluid dynamics. Magnon wave scattering off a localized Skyrmion is examined numerically based 
on our Gross-Pitaevskii formulation. Dependence of the effective flux experienced by the impinging 
magnon on the Skyrmion radius is pointed out, and compared with analysis of the same problem 
using the Landau-Lifshitz-Gilbert equation. 

PACS numbers: 75.78.n, 75.30.Ds, 12.39.Dc, 75.78.Cd 


I. INTRODUCTION 

Magnons are quantized waves of spin fluctuations 
around an ordered state of spins in a magnet. The theory 
of their dynamics was formulated in the 1940s by Holstein 
and Primakoff [1] and for this reason they are also known 
as Holstein-Primakoff (HP) bosons. Magnon dispersion 
in the ferromagnetic or the anti-ferromagnetic spin back¬ 
ground can be worked out easily from the HP theory as 
explained in textbooks on quantum magnetism [2]. 

Modern experimental and theoretical advances call 
for refinements of the HP theory. On the experimen¬ 
tal front abundant sightings of the Skyrmionic spin tex¬ 
ture in chiral magnets demand a well-defined theory of 
spin waves in the textured spin background [3]. Exist¬ 
ing attempts [4-8] for theories of spin dynamics in the 
Skyrmion background are done in terms of the Landau- 
Lifshitz-Gilbert (LLG) equation for the magnetization 
unit vector n(r, t). Despite its obvious strength in direct 
simulation of the magnetization dynamics in such com¬ 
plex spin background, LLG approach is not well suited 
for an intuitive understanding of the low-energy spin dy¬ 
namics based on the particle picture. Moreover, magnon 
scattering off the localized Skyrmion is a subject of grow¬ 
ing importance and interest as addressed by several re¬ 
cent works [7, 8]. Here again, a formulation of the prob¬ 
lem in terms of the well-established HP boson theory is 
surprisingly lacking. 

As amply demonstrated in several recent theories of 
magnon Hall effect [9-13], HP theory conveniently cap¬ 
tures magnon dynamics in a non-trivial spin background 
in a manner completely analogous to that of electron dy¬ 
namics in a non-trivial flux background and band Chern 
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number. In these approaches, however, the relevant de¬ 
scription is within the momentum space picture assuming 
translationally invariant spin background over some large 
unit cell. We try in this article to present a thorough, self- 
contained formulation of magnon dynamics in real space , 
suited in particular for adressing the scattering process 
off a localized object such as a Skyrmion. The work is 
thus complementary to both strands of existing theories 
of spin dynamics based on LLG equation [7, 8], and those 
formulated in momentum space [9-13]. 

In section II we present a general HP Hamilto¬ 
nian assuming arbitrary smooth spin background of 
the ferromagnetic exchange Hamiltonian supplemented 
by Dzyaloshinskii-Moriya (DM) exchange and Zeeman 
terms. While this formulation was attempted in some 
earlier works [14, 15] and applied in various context of 
spintronics, not all of the terms we derive here have been 
coherently presented thus far. In particular the terms di¬ 
rectly following from the DM spin interaction have never 
been consistently derived to our knowledge. Then in sec. 
Ill we show how the magnon equation of motion pre¬ 
sented in the previous section can be equivalently viewed 
as a hydrodynamics problem, with a set of hydrody¬ 
namic equations pertaining to magnon density and veloc¬ 
ity flows. Numerical simulations based on HP theory are 
in fact equivalent to solving the hydrodynamics version of 
the magnon problem. We present such numerical result 
in sec. IV, in the particular case of magnon scattering off 
a single, localized Skyrmion. Our results partially over¬ 
lap with an earlier simulation of the same problem [8]. 
Here we find a surprising dependence of the magnon scat¬ 
tering angle (and even its sign!) on the Skyrmion radius 
normalized by the spiral wavelength. This unexpected 
feature of the magnon-Skyrmion scattering follows from 
taking thorough consideration of the DM term in the HP 
formulation. We carry out a numerical check of the pre¬ 
diction using the more conventional LLG approach. It 
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turns out that LLG approach embodies a more complex 
picture of the scattering dynamics due to the fact that the 
Skyrmion tends to respond dynamically to the impinging 
magnon wave in the LLG theory, whereas the HP theory 
in effect treats it as a static object of infinite mass and 
internal excitation energies. Despite these differences, we 
were able to deduce features of the subtle dependence of 
the magnon scattering on the Skyrmion radius that was 
unexpected in previous simulations [7, 8]. We conclude 
with a summary and possible future applications of our 
formulation in sec. V. 


where is the unit vector in each /r-direction. Hence¬ 
forth we choose the unit of energy J — 1 and re-scale 
the field strength B accordingly. Restriction to d = 2 
is easily achieved by deleting all terms from Eq. (2.4) 
pertaining to the spatial direction z. 

Let us denote the new ground state or some metastable 
spin configuration of the Hamiltonian H by no- Given 
such 


no = (sin 0 cos </>, sin 0 sin 0, cos 0) 
one can construct an orthogonal matrix R whose elements 
are 


II. HAMILTONIAN FORMULATION OF 
MAGNON DYNAMICS 

The low-energy spin dynamics of a ferromagnet is cap¬ 
tured by the nonlinear sigma model (NLcrM) 


Ho = \ ■ (0„n) (2.1) 

fl=l 

with the spatial index /i running over 1 through d in 
d -dimensional space. When a smooth perturbation is 
added to the NLcrM the new ground state will gener¬ 
ally become a slowly varying spin texture. In light of the 
recent surge of interests in textured ground states of the 
chiral magnet [3], we include the DM interaction, 


Ha/3 — 2777 / q , TYlf3 da(3i 

, Q Q Q v 

m = ^ sin - cos 0, sin - sin 0, cos - j . (2.5) 

With this R one can show Rz = no and Rno = z. The 
strategy is to introduce HP bosons in the rotated frame 
of reference, as obtained by replacing n in Eq. (2.4) by 

Rn. 

After some extensive algebra one finds the Hamiltonian 
in the new spin basis 


l u 

H = - ^2 -a M xnf-SBn 0 -n (2.6) 
1 

where, component-wise, the vector potentials are 


Edm = D n • V x n, (2.2) 


as well as the Zeeman term 


Hz = -SB • n (2.3) 

(.S = size of spin) as the perturbation. Space integration 
is implicitly assumed in the above expressions. It is useful 
for further technical development to point out that the 
combined spin Hamiltonian H = Ho + Hbm + Hz can 
be organized in the compact form, up to an irrelevant 
constant and taking k = D/J, 


H 

7 


~ E - Ki n x n ) 2 

^=i 


s 

J' 


Bn, 


(2.4) 


a l = ~ n x ( < d ^ ( t ) ) - sin<£(d M 0) + k/ m (6»,(/>), 

al = cos (j){d^9) - riyid^) + ng^{e,(f>), 

a l = (! - cos9)(d ll (j)) + kh“. (2.7) 


The two functions ^(0,0) and ^(0,0) tied to the DM 
interaction are 




9 ^( 6 , 4 > ) 


— cos 2 | + cos 2 (j) sin 2 | 

sin 2 cj) sin 2 | 

sin 2 (j) sin 2 | 

— cos 2 | — cos 2 <p sin 2 | 


(g = x) 

(m = y) > 

(m = z) 

(m = x) 

b = y) ■ ( 2 - 8 ) 

(M = z) 


All the /^-dependent terms in the emergent gauge fields 
were missed in the literature. 

Quadratic magnon Hamiltonian follows from the HP 
substitution in the rotated spin Hamiltonian (2.6): 
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H 

~S 


d d 

o y [(^+*<4)4] [( d K~ ia3 M + 9 J2i^~ ia3 M [( d »+ ia D b l\ 


ii=i 


m=i 


B • n 0 - \ ^[( a M) 2 + ( a i^) 2 ]) b r b r ~ 7 ^2( a l~ ia l) 2b r b v ~ T E] (a^ + M^) 2 ^! 


M=1 


M=1 


M=1 


(2.9) 


The size of spin S that serves as an overall constant in the 

_I 


quadratic theory can be set to one. Heisenberg equation 
of motion for the boson operator is 


<%r 

lh ~iH = ~ - « 

11=1 


3 \2 


)% 



This equation of motion for the magnon operator 6 r , com¬ 
bined with Eqs. (2.7) and (2.8) for the background gauge 
fields, gives the governing dynamics of magnons moving 
in the arbitrary smooth spin background no. Relation of 
the gauge field a 3 to the background spin n 0 is 

(V x a 3 ) a = no • (d M n 0 x d v n 0 ) + k(V x n 0 ) Q 
+(1 — cos0)(V x V0) a . (2.11) 

The prefactor in the first term on the r.h.s. shows that 
one Skyrmion acts as two units of flux quanta for the 
bosons. Interestingly, the spin vorticity V x no also con¬ 
tributes to V x a 3 due to the DM interaction. We will 
see later that this extra contribution has an interesting 
consequence for the magnon-Skyrmion scattering prob¬ 
lem. The final expression on the r.h.s. might be ignored 
on the ground that V x V</> = 0, but as we will soon 
see this is not the case with the Skyrmion or the anti- 
Skyrmion configuration we have in mind. 

The anomalous potential, which violates the boson 
number conservation, has the connection to the back¬ 
ground spin 


- Y(d+*4) 2 


ii=i 


= e 2 ^ 


a 2 

y (^+^)-(^n 0 -Ke M xno) -(2-12) 


11=1 


Other unit vector fields in the above formula are 


6 = (cos 0 cos 0, cos 0 sin </>, — sin 0 ), 
(j) = (— sin 0, cos 0, 0), 

both locally orthogonal to n 0 . 



[K? + (a') 2 ] 


( 2 . 10 ) 


The meaning of the anomalous term becomes clear 
with the exemplary spiral spin no = (cos ft 2 , sinftz, 0) 
- a right-handed spiral with the propagation vector k = 
(0,0, ft). We find that Eq. (2.12) reduces to ft 2 e 2m2 , and 
the whole magnon Hamiltonian written in momentum 
space becomes 



One recovers the well-known magnon dispersion of the 
spiral spin, cj q = |q|y/q 2 + ft 2 . It is through the anoma¬ 
lous potential that the correct magnon dispersion is re¬ 
covered, while its omission might have led to an incorrect 
dispersion c<; q = q 2 + ft 2 /2. 

Finally, there is an inhomogeneous on-site potential 

y [(fl£) 2 + (a-l) 2 ] = y {df.no -x n 0 ) 2 (2.14) 
11=1 / 1=1 

arising from local deviation of the spin structure from 
that of a simple spiral. 

The formulation presented in this section is com¬ 
pletely general, and applies to arbitrary ground state or 
metastable spin configurations allowed in the Hamilto¬ 
nian (2.4). 


III. HYDRODYNAMIC FORMULATION OF 
MAGNON DYNAMICS 

Magnon dynamics is typically viewed as the solution 
of Eq. (2.10) with some characteristic frequency cor¬ 
responding to the energy of magnon excitation. In an 
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inhomogeneous environment of the spin no, however, it 
will be useful to have a complementary picture in real 
space, quite like that provided by the hydrodynamic for¬ 
mulation of superfluid dynamics [16], which follows from 
writing down the corresponding magnon Lagrangian and 
substituting 

b = (3.1) 

Variation of the Lagrangian with respect to the density 
p and the phase p yields hydrodynamic equations. An¬ 
other way, paralleling the derivation of Gross-Pitaevskii 
equation in superfluids [16], is to replace the boson op¬ 
erator b by its average ( b ) = y/pe ir] in the Heisenberg 
equation of motion (2.10). Same hydrodynamic relations 
are obtained from both approaches. 

The boson number is no longer conserved in the con¬ 
tinuity equation due to the anomalous potential, 

^p + V ■ (pv) = pFi, (3.2) 

where Fi follows from 


the real-space magnon dynamics, presumably based on 
the observation that they appear in the magnon Hamil¬ 
tonian with one more spatial derivative than a 3 . On the 
other hand, we already saw that their omission leads to 
an erroneous result even in the well-known test case. In¬ 
stead the relevance of the anomalous term has to be ad¬ 
dressed carefully, for each specific background spin con¬ 
figuration. 

In the next section we take up in particular the prob¬ 
lem of magnon scattering off a single localized Skyrmion 
based on our HP formulation. This problem is of great 
topical interest and has already received some attention 
in the literature [7, 8]. The direct numerical approach 
based on LLG equation in these papers however makes 
it difficult to interpret the scattering process in terms of 
magnon particles. Solving the scattering problem analyt¬ 
ically together with the anomalous term is impossible 8. 
Our method of choice is the numerical integration of the 
magnon dynamics, as given in Eq. (2.10). 

IV. SIMULATION OF MAGNON SCATTERING 
A. Simulation based on HP theory 


~ e ~ 2iV (xx+^) 2 J = p r + iFi- (3-3) 

The other piece of hydrodynamics is provided by the Eu¬ 
ler equation of the flow vector v = 2(Vr? — a 3 ) given 
by 

D t w = vx8 - Vp. (3.4) 


For numerical simulation based on Eq. (2.10) we work 
in two dimensions, external field B = Bz , B > 0, and 
the DM interaction constant k > 0 that favors a right- 
handed spiral. Under these parameter values the stable 
Skyrmion configuration is the one whose spins point up 
far from the core, execute a right-handed spiral turn upon 
approach to the origin, where it points downward, —z. A 
model spin configuration fulfilling all these requirements 
is (r 2 = x 2 + y 2 ) 


The material derivative D t = dt + v • V is familiar from 
hydrodynamics. The emergent magnetic field B respon¬ 
sible for the Lorentz force on the magnon is defined by 

B = —^ V xv = Vxa 3 — Vx Vp. (3.5) 


2 Ry 2 Rx r 2 - R 2 

r 2 + R 2 ' r 2 +R 2 ' r 2 + R 2 

The Skyrmion number of this configuration is N = — 1, 
qualifying it as an anti-Skyrmion. The HP Hamiltonian 
in the presence of a single anti-Skyrmion is 




Readers can refer to Eq. (2.11) for the definition of V x 
a 3 . We will not immediately identify B with the curl 
V x a 3 because in some situations the phase singularity 
in the boson field can lead to V x Vp ^ 0. Allowing for 
the time-dependent background spin (which we do not 
do in this paper) will induce the emergent electric field 
on the r.h.s. of the Euler equation as well. The quantum 
pressure p is given by 


P = ~ 


—-=-^ + 2 BriQ + F r -^2 [«) 2 + (^) 2 ] -( 3 - b ) 

VP n=i 


2 

H = + «4)4][(<^ - ia l) h r] - G{r)b\b r 

M=1 

+l e - ii +F[r)b T b T + le 2< *F(r)&J&J, (4.2) 

where 


F(r) = 


G(r) = 


2 kR(kR - 2)r 2 
(. R 2 + r 2 ) 2 

R 2 (kR — 2) 2 + ft 2 r 4 R 2 — r 2 
(R 2 +r 2 ) 2 R 2 + r 2 


(4.3) 


There is a tendency in the existing literature [8, 14, 15, 
17, 18] to ignore effects arising from a 1 and a 2 in viewing 


The a 3 and its curl is tricky as it contains a singular 
contribution, 
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3 _ ( -2 yR 2 2xR 2 \ 

a y r 2( r 2 _|_ ’ r 2(^ r 2 J 

f —2 Ry 2 Rx \ 

\r 2 + R 2 ’ r 2 + R? / 

( V X a ), = (r 2 + p>2)2 + 4 ^W- ( 4 - 4 ) 

The delta-function flux can be removed through the sin¬ 
gular phase transformation of the boson field b r —> b r e 
where <j) is the azimuthal angle in the plane. The new bo¬ 
son field obeys the Hamiltonian 

2 

H = + ia l) b t][( d n - iafybr] - G(r)blb r 

n=l 

+ ^e 2i *F(r)b r b T + le~ 2i *F(r)bibl (4.5) 

where F and G remain the same as before, but a 3 and 
its curl becomes 


a 


3 


(V x a 3 ) z 


( 2y( K R-l) 
\ r 2 + R 2 
AR 2 (kR — 1) 
(r 2 + R 2 ) 2 


2x(kR-1)\ 
r 2 + R 2 J 


(4.6) 


In the numerical treatment we will work in this 
singularity-free gauge. 

It is customary to regard the two-dimensional integral 
of (V x a 3 )^ divided by 2tt as the effective flux experi¬ 
enced by the magnon quasiparticle. In such picture each 
Skyrmion serves as 


$ = 2(kR — 1) 

units of flux quanta, spread over the distance of the ra¬ 
dius R. In contrast to electron motion which sees the 
Skyrmion as a quantized flux object carrying one flux 
quantum [3], magnons see it carrying a variable flux that 
can even change sign with the Skyrmion radius! 

Let us rescale all the fields into dimensionless form 
convenient for numerical analysis, 


F(f) = 


G(f) = 


a 3 = 


2 k(k — 2 )r 2 
(1 + r 2 ) 2 ’ 

(k - 2) 2 + K 2 f 4 - 1-f 2 

(1 + f 2 ) 2 + 1 -hr 2 ’ 

/ 2 y(k - 1) 2 x(k - 1) \ 

\ 1 + f 2 ’ 1 + r 2 / ’ 


(4.7) 


where x = x/R, y = y/R, r = r/R , and k = kR and 
B = BR 2 are the two dimensionless parameters. The 
equation of motion in dimensionless form becomes 


(a) (b) (c) 



FIG. 1. (color online) Snapshots of scattered magnon waves 
for k = kR equal to (a) 0.5, (b) 1.0, and (c) 1.5. Shown are 
color plots of the magnon density p — |6| 2 normalized by the 
incoming magnon density po. Skyrmion is indicated by a red 
circle of radius R. The magnon trajectory changes direction 
at k = 1 as anticipated from the total flux consideration in 
the HP theory. 



(<9 m - fa 3 ) 2 + G{f) 


h + e~ 2i<l, F(f)b*„ 


(4.8) 


where the rescaled time t = t/h is also dimensionless 
(Recall that we already set J = 1 earlier). The tilde can 
be removed now without confusion. 

Hydrodynamic relations following from Eq. (4.8) are 


+ V- (pv) = —2pF sin[2(?7+0)], (4.9) 

and 

D t w = v x B — Vp, 

p = 2 ^-- G + Fcos [ 2 ( 77 + 0 )]^ . (4.10) 

We solved Eq. (4.8) using the Runge-Kutta method. 
Forced spin wave is generated at the bottom of the rect¬ 
angular simulation grid with the space-time dependence 


b(x,y<y 0 ,t)=b 0 e ik y- i ^ t . (4.11) 

The frequency c Ok = k 2 +B matches the dispersion in free 
space far from the Skyrmion, and yo is a suitable cutoff 
much less than the total length of the grid. We then 
closely monitor the scattering processes as the magnons 
collide with the (static) Skyrmion. 

By far the most interesting findings of our investiga¬ 
tion was the variation of the scattering angle as k = kR 
crosses the threshold value 1, whereupon the scattering 
direction changes sign as shown in Fig. 1. This is a stark 
prediction of our HP theory, in departure from previous 
considerations based either on LLG theory [7, 8] or on 
the magnon theory without taking into careful account 
of the DM term [14, 15]. In the following subsection we 
return to the conventional LLG calculation to verify to 
what extent, if at all, the predictions of the HP theory is 
born out. 
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B. Simulation based on LLG theory 

Lattice version of the LLG equation, 

hi = UiX Bf - an * x (n* x B- ff ) , 

B® ff = ^ (rij + ctKUj x ej) + a 2 B, (4.12) 
jei 

is solved on a large grid of lattice spacing a, with the 
initial spin configuration given by Eq. (4.1). It was found 
useful to first run the simulation with a large Gilbert 
constant a of order one to relax the Skyrmion to its true 
equilibrium configuration, which differs in detail from the 
model form (4.1), at a given B and k. Once an optimal 
Skyrmion configuration is reached in this way, we turn 
on the spin wave at the bottom of the boundary to start 
the scattering process. 

Several different Skyrmion radii R were considered 
in the simulation while simultaneously maintaining the 
same ratio to the magnon wavelength, kR. A detailed de¬ 
pendence of the spin wave scattering angle on this dimen¬ 
sionless quantity kR was reported earlier [8], whereas our 
focus here is with another another dimensionless number, 
kR. Variation of this number can be achieved in practice 
with the magnetic field R, which would choose a different 
optimal Skyrmion radius R. 


(a) (b) (c) 



(d) (e) (£) 



FIG. 2. (color online) Snapshots of scattered spin waves in 
the LLG calculation for kR equal to (a,d) 0.9, (b,e) 1.5, and 
(c,f) 2.0. In the top row are plots of the spin wave den¬ 
sity p — yj(6n x ) 2 + ( Sn y ) 2 . See the main text for the def¬ 
inition of S n = ( Sn x ,Sn y ). Bottom row shows the phase 
4> = arctan($n y /#n x ). Skyrmion positions and their sizes are 
indicated by red circles. A clear sign of skew scattering at 
kR = 0.9 is completely gone at kR = 2.0. 

Figure 2 shows patterns of scattered spin waves for 
varying Skyrmion radii R. The background spin texture 
has been removed by taking the time average (rq) over 
a single period of the incoming spin wave, which is then 
subtracted from the snapshot at a particular instant rq 


to obtain $rq = rq — (rq). The resulting profile has only 
planar components <5rq = (. 5n X: Sn y ) whose magnitude 
and phase are plotted in Fig. 2. A small amount of 
displacement of the Skyrmion [8] during the cycle does 
not affect the average (rq) significantly. 

At the smaller radius kR = 0.9, overall scattering be¬ 
havior are analogous to that of magnons for kR < 1. 
At kR = 1.5, we see a second scattering channel begin¬ 
ning to open up shown as a white, vertical jet in Fig. 
2(b). At kR = 2.0, this second channel is pointing at 
an angle opposite to the first one (Fig. 2(c)), while the 
phase profile (Fig. 2(f)) becomes nearly isotropic with 
respect to the incoming spin wave direction. An attempt 
to carry out LLG calculation at even larger radius R, 
in the hope of finding the reversed scattering direction, 
failed because stabilizing such a large Skyrmion requires 
B for which Skyrmion is no longer the most stable phase 
(instead spiral spin state wins out). To further compli¬ 
cate the analysis, low-energy excitation modes are easily 
created for larger Skyrmions which subsequently interact 
with the incoming spin wave. 

In short, our LLG simulation does give a hint that 
spin wave scattering angle has a significant dependence 
on the reduced Skyrmion radius kR, in accordance with 
anticipations of the HP theory. To observe the proposed 
effect one can tune the magnetic field over a range in 
which isolated Skyrmions are the most favored states [3]. 
Field-dependent variations of the magnon Hall effect and 
possible reversal of its sign will be a telltale sign of our 
effect. 


V. SUMMARY AND DISCUSSION 

Much of this paper dealt with the formal theory of 
magnon dynamics from Hamiltonian and hydrodynam¬ 
ics perspectives. Despite the long history of the theory 
of magnon dynamics dating back to 1940s, its real-space 
formulation with regard to the underlying spin Hamil¬ 
tonian containing the Dzyaloshinskii-Moriya interaction 
has never been rigorously derived in the form presented 
in Sec. II. 

The spin Hamiltonian (2.4) used in this paper is known 
to host interesting topological phases of spins called 
the Skyrmion lattice in two dimensions [19, 20] and the 
hedgehog lattice in three [21, 22]. The dynamical equa¬ 
tion derived in this paper can be used for large-scale sim¬ 
ulation of magnon dynamics in these topological phases. 
The hydrodynamic formulation developed in Sec. Ill can 
subsequently provide intuitive interpretation of magnon 
dynamics simulation as movements of magnon fluid. 

In this paper we studied one such problem - magnon 
scattering off a localized Skyrmion - numerically in de¬ 
tail. Magnons do see a localized Skyrmion as a source 
of flux quanta responsible for skew scattering, in accor¬ 
dance with several earlier observations [7, 8, 15], but the 
net flux seen by the magnon deviates from the two units 
of flux quanta anticipated in earlier theories due to an 
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important correction from the DM interaction, as in Eq. 
(2.11). Our analysis suggests that the net flux seen by 
the incident magnon may even change its sign at the crit¬ 
ical Skyrmion radius kR c = 1. Indeed magnon dynam¬ 
ics simulation based on our HP Hamiltonian did revealed 
such effect. Unfortunately, the LLG simulation presented 
a more complex picture due to various low-energy exci¬ 
tations of the Skyrmion induced by the very magnons 
which scatter off of it. Nevertheless we were able to 
show that a Skyrmion of increasingly larger radius has 
less pronounced skew scattering of magnons in rough ac¬ 
cordance with the prediction of the HP theory. We antic¬ 
ipate that the magnon Hall effect due to Skyrmions will 
have a more complex dependence on the external mag¬ 
netic field (which controls the Skyrmion radius to some 
extent) than the topological Hall effect experienced by 


the electrons. 
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